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ABSTRACT 


Seismic  event  locations  obtained  by  using  P  and  P  arrivals  may  cause 

g  n 

severe  errors  because  the  Earth's  crust  is  laterally  heterogeneous.  Location 
errors  are  attributed  to  the  fact  that  crustal  velocities  and  thicknesses  for 
each  source-to-station  path  are  different  from  those  of  the  standard  crustal 
model  used  in  the  location  program.  To  overcome  this  difficulty,  we  have 
installed  a  number  of  crustal  models  in  a  location  program  and  made  it  possible 
for  each  station  to  be  assigned  an  appropriate  model.  Travel  times  for  each 
phase  are  computed  by  using  the  selected  earth  model  and  the  specified  phase. 

The  merit  of  this  method  is  evaluated  by  relocating  12  explostions  in  the 
western  United  States.  When  each  station  model  was  chosen  individually  from 
a  number  of  regional  model.,,  location  errors  were  slightly  larger  than  errors 
using  only  one  (regional)  model  for  all  stations.  However,  errors  were  reduced 
when  each  station  was  assigned  its  ow  .  model  from  a  set  of  finely  localized 
crustal  models. 

Errors  were  also  reduced  when  the  local  crustal  model  appropriate  to  the 
source  region  was  used  for  all  stations.  This  suggests  that  a  crustal  model 
for  the  source  and  a  separate  model  for  each  station  would  result  in  even 
better  locations. 

The  addition  of  P  arrivals  to  the  set  of  P  arrivals  reduced  location 
g  n 

error  by  about  30a.  However,  the  difference  is  not  great  enough  for  this  small 
sample  but  we  can  be  sure  that  it  is  significantly  different  from  zero. 
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1.  INTRODUCTION 

The  accuracy  of  event  location  depends  strongly  upon  the  variation  of 
crustal  structures  from  the  seismic  source  to  each  station.  We  have  devised 
a  method  to  minimize  the  effect  of  crustal  variations  by  using  various  crustal 
models  for  each  source-to-station  path. 

Chiburis  and  Ahner  (1970)  studied  the  location  accuracy  of  28  Nevada  Test 
Site  explosions  using  teleseismic  P  arrivals.  They  showed  that  the  location 
accuracy  of  these  events  (depth  restrained)  averaged  about  7.66  kilometers 
but  could  be  as  large  as  20  kilometers.  They  repeated  the  relocation  of  those 
events  with  station  corrections  and  location  errors  were  reduced  to  about 
2.81  kilometers.  However,  Chiburis  and  Ahner  cautioned  that  the  set  of  station 
corrections  for  the  Nevada  Test  Site  was  applicable  only  to  a  small  region. 

Location  accuracies  resulting  from  using  P  arrivals  at  local  stations 
were  investigated  by  Engdahl  and  Lee  (1976).  Three  methods  were  compared  in 
that  study.  The  first  method  was  HYP071  (Lee  and  Lahr,  1972),  where  a  uniformly 
layered  crustal  model  was  used.  The  second  method  was  HYP074,  where  different 
crustal  models,  as  well  as  station  corrections,  were  used  for  each  individual 
station.  The  third  method  used  a  complex  two-dimensional  model  to  describe 
crustal  variations  across  the  San  Andreas  fault  near  Bear  Valley.  Travel  times 
for  each  station  in  the  third  method  were  computed  by  a  ray  tracing  method. 
Although  the  third  method  is  the  most  elaborate  one,  improvements  in  location 
accuracies  over  the  second  method  appeared  to  be  small.  All  stations  used  in 
their  study  were  within  30  kilometers  of  the  epicenter. 


Chiburis,  E.  F.  and  R.  0.  Ahner  (1970).  A  seismic  location  study  of  station 
anomalies,  network  effects,  and  regional  bias  at  the  Nevada  Test  Site, 
Teledyne  Geotech,  Seismic  Data  Laboratory  Report  No.  253,  Alexandria,  VA. 

Engdahl,  E.  R.  and  W.  H.  K.  Lee  (1976).  Relocation  of  local  earthquakes  by 
seismic  ray  tracing,  J.  Geophys.  Res.,  81,  4400-4406. 

Lee,  W.  H.  K.  and  J.  C.  Lahr  (1975).  HYP071  (Revised):  A  computer  program 
for  determining  hypocenter,  magnitude,  and  first  motion  pattern  of 
local  earthquakes,  USGS  Open-file  report  75-311,  Menlo  Park,  CA. 
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2.  METHOD  OF  INVESTIGATION 


2.1  The  Location  Program 


The  location  program  used  in  this  project  was  originally  developed  by 

Julian  (1974).  McCowan  (1978)  modified  it  and  installed  an  option  to  compute 

travel  times  for  six  crustal  phases,  P  ,  P*,  P  ,  S  ,  S*,  and  S  ,  from  a  crustal 

n  g  n  g 

model  ( Jef f reys-Bullen  model) .  We  modified  the  program  to  allow  it  to  use  up  to 
fifty  crustal  velocity  models.  Furthermore,  the  program  was  modified  such  that 
each  station  can  use  an  independent  model  to  compute  travel  times. 


Crustal  Models  and  Their  Usage 


Figure  1  shows  the  boundaries  and  codes  for  ten  regional  crustal  models 
published  by  Pakiser  and  Robinson  (1966),  of  which  numbers  1  through  7  are  used 
in  this  investigation.  In  Figure  2,  fifteen  more  localized  crustal  models  in 
the  southwestern  United  S cates  (SWUS)  are  shown.  In  addition  to  those  models, 
crustal  models  for  Herrin  68  and  Jef freys-Bullen  are  also  available  in  the 
program.  Table  I  is  a  comparative  listing  of  these  models. 


The  input  to  the  location  program  must  identify  the  station,  signal  phase, 
and  the  crustal  model  to  be  used  for  this  station.  To  identify  the  crustal 
model  the  analyst  must  make  a  preliminary  location  of  the  source  and, 
with  Figures  1  and  2,  decide  for  each  station  on  the  best  overall  model  for  the 
path  from  the  source  to  the  station.  Although  this  scheme  provides  versatility, 
it  is  sometimes  difficult  to  choose  a  model,  as  a  particular  source-to-receiver 
station  path  may  travel  across  regions  best  described  by  several  different 
models.  In  such  a  case  the  analyst  may  opt  to  choose  an  average  model,  such  as 
the  Herrin  model,  for  this  station. 


Julian,  B.  (1973).  Extension  of  standard  event  location  procedures.  Seismic 
Discrimination  SATS,  Lincoln  Laboratory,  M.I.T.,  30  June  1978,  4-9. 

McCowan,  D.  W.  (1978).  Personal  communication. 

Pakiser,  L.  C.  and  R.  Robinson  (1966).  Composition  of  the  continental  crust 
as  estimated  from  seismic  observations.  The  Earth  Beneath  the  Continents. 
American  Geophysical  Union,  Geophysical  Monograph  #  10,  620-626. 
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TABLE  I 


Comparison  of  Crustal  Structures 


1st  Layer 


2nd  Layer 


Mantle 


Area  I  Designator  Thickness  Velocity  Thickness  I  Velocity  Velocity 


General: 

Jeffreys-B ullen 
Herrin  68 

Regional; 

Calif.  Coast 
Sierra  Nevada 
Pac.  NW  Coast 
Columbia  Plat. 
Basin  &  Range 
Colorado  Plat. 
Rocky  Mtns. 

Local: 

N.  Calif. 

Coast  Calif. 
Sierra  Nevada 
S.  Calif. 

N.  Nevada 
Cent.  Nevada 
SW  Nevada 
Lake  Mead  Nev. 
W.  Utah 
E.  Utah 
N.  Arizona 
Cent.  Arizona 
S.  Arizona 
W.  Colorado 
W.  New  Mexico 


4 


2.3  Data 


Twelve  nuclear  explosions  were  used  as  the  data  base.  Of  these  events 
seven  were  at  the  Nevada  Test  Site  (NTS)  and  six  others  were  located  in 
Colorado,  New  Mexico  and  Nevada. 


and  P^  arrivals  for  four  events  were  picked  by  analysts  at  the  Seismic 
Data  Analysis  Center  with  cross  checks  to  the  published  shot  reports  for  each 
event.  For  the  rest  of  the  events,  phase  arrivals  as  reported  in  the  shot 
reports  were  used.  Reading  errors  of  and  P^  are  greater  than  errors  for 
teleseismic  P.  P^  and  P^  (especially  P  )  signal  waveforms  are  very  complex 
and  often  emergent.  In  Table  II  we  list  the  event  parameters,  numbers  of  P^ 

and  P  arrivals,  and  the  data  sources  for  the  events  used. 

g 

2 . 4  List  of  Location  Experiments 


In  the  following  test  categories  events  were  relocated  using  (a)  P^  and 
P  phases,  (b)  P^  only,  and  (c)  P^  only: 

1)  Local  station  models; 

2)  Regional  station  models; 

3)  Herrin  68; 

4)  Epicenter  models;  and 

5)  Jef freys-Bullen. 

The  Jef freys-Bullen  model  was  not  tested  with  P^  only,  because  large  re¬ 
sidual  times  were  observed  using  this  model.  Location  errors  using  P^  only 
were  so  large  that  it  was  eliminated  from  further  analysis.  It  was  necessary 
to  restrict  the  event  depth  to  zero  in  all  tests,  because  the  hypocenter  depth 
may  go  deeper  than  the  crustal  thickness  during  iterations  in  which  case 
P  and  P  travel  times  can  not  be  computed. 

n  8 
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3.  DISCUSSION 


In  Table  III  location  errors  in  latitude  and  longitude  are  combined  to 

give  errors  in  kilometers  from  the  true  location.  The  entries  in  the  column 

giving  the  best  model  are  determined  by  consideration  of  location  error,  origin 

time  error,  and  station  residuals.  Origin  time  errors  in  the  form  of  T  ,  -T 

calc  true 

are  given  in  Table  IV. 

3.1  Location  Accuracy  with  P  and  P 
_ 2 _ g. 

Comparing  location  and  origin  time  errors  of  runs  with  P^  +  P  and  with 
P^  only,  one  notes  that  the  location  accuracy  with  Pg  is  only  about  30%  better. 
The  result  is  not  clearly  significant  and  suggests  that  the  addition  of  to  the 
existing  P^  arrivals  may  not  help  to  improve  location. 

The  best  result,  an  error  of  6.31  km,  was  obtained  using  crustal  epi¬ 
center  models  and  ?n  +  P  •  The  average  origin  time  error  1.31  seconds  is 

high  for  this  model.  The  average  errors  using  the  Herrin  model  and  P  arrivals 

n 

are  8.41  km  and  0.53  sec.  The  average  location  errors  in  Table  III  show  that, 
using  proper  models  the  accuracy  of  event  location  with  P^  and  P^  is  about  6 
to  8  kilometers,  or  approximately  equal  to  the  location  accuracies  using  tele- 
seismic  P. 

Average  location  errors  using  regional  models  (10  to  11  km)  are  higher 
than  errors  using  an  epicenter  model  (6  to  8  km).  Location  errors  using  local 
models  (7  to  8  km)  are  comparable  to  average  errors  using  the  Herrin  model 
(-  8  km)  but  with  better  origin  time  errors. 

These  results  show  that  an  accurate  source  model  which  is  correct  for 
portions  of  all  paths,  or  local  station  models  which  correctly  handle  the 
upcoming  ray  at  each  station,  give  good  results.  The  implication  is  that 
several  models  are  needed  for  each  path  to  attain  excellent  results:  a 
source  model  for  the  down-going  ray,  a  receiver  model  for  the  up-coming  ray, 
and  a  path  model  to  give  the  average  propagation  velocity.  However,  the 
present  program  is  designed  for  only  one  model  per  path. 

In  an  attempt  to  minimize  the  effect  of  crustal  variations,  Herrin  and 
Taggart  (1962)  computed  an  average  P^  velocity  and  an  average  crustal  thick¬ 
ness  for  each  individual  path.  The  crustal  velocity  was  assumed  fixed  through- 


TABLE  III 


applicable  because 


TABLE  IV 


applicable  because  local  models  not  available  for  all  stations 


out  the  path.  We  feel  that  this  does  not  yield  a  complete  solution  to  the 

problem  because  travel  time  variations  due  to  crustal  thickness  at  source  and 

receiver  are  much  greater  than  travel  time  variations  due  to  differences  in 

velocity.  For  example,  assume  a  particular  source-to-receiver  path  across 

the  adjacent  Basin  and  Range  (B  &  R)  and  Colorado  Plateau  (CP)  regions  where 

the  ray  travels  through  250  km  of  B  &  R  and  250  km  of  CP  upper  mantle.  The 

difference  in  P  velocities  between  these  regions  is  .1  km/sec,  i.e.  7.9  for 
n 

B  &  R  and  7.8  for  CP  (Table  I);  therefore,  if  one  model  were  used,  the  error 
in  P^  travel  time  would  be  0.4  second.  On  the  other  hand,  using  the  thick¬ 
nesses  of  the  B  &  R  crust  (30  km)  and  the  CP  crust  (40  km)  and  an  average 
crustal  velocity  of  6.6  km/sec  for  both  (Table  I),  it  can  be  seen  that  the 
shortest  one-way  surface-to-mantle  travel  time  is  1.6  seconds  longer  for  the 
CP  than  for  the  B  &  R.  Thus  using  one  model  in  this  case  would  mean  an  error 
of  at  least  1.6  seconds  in  the  crustal  travel  time  while  the  error  resulting 
from  the  incorrect  P^  velocity  would  be  a  quarter  of  that  at  0.4  second. 

This  example  demonstrates  that  the  effect  of  the  P^  velocity  on  location 
accuracy  is  smaller  than  the  effect  of  the  variations  in  crustal  structures. 
This  also  demonstrates  why  using  regional  models  appropriate  to  the  path  but 
not  to  the  crustal  structure  under  the  station  may  not  improve  locations. 

3 . 2  Summary  and  Recommendations  for  Future  Research 

The  location  method  used  in  this  study  was  only  marginally  successful 
in  improving  location  accuracy.  This  is  mainly  due  to  the  heterogeneity  of 
the  crust  in  the  western  United  States  (WUS).  The  method  should  produce 
better  results  in  an  area  where  lateral  crustal  variations  are  not  as  severe 
as  those  in  the  WUS.  The  reported  success  in  calibrating  teleseismic  P  with 
the  master  event  method  (Chiburis  and  Ahner,  1970)  may  not  hold  with  respect 
to  locations  with  crustal  phases  because  of  the  stronger  effect  of  lateral 
heterogeneity  on  these  phases.  How  accurately  the  models  we  used  reflect 
actual  crustal  conditions  is  open  to  question  but  this  is  a  matter  outside 
the  scope  of  this  study. 

P^  observations  proved  not  to  be  as  helpful  as  expected  in  improving 
location  accuracy  perhaps  because  the  onset  of  this  phase  is  obscured  by  coda. 
Consequently,  and  P^  should  not  be  weighted  equally  in  computing  locations 
but  each  observation  should  be  weighted  as  a  function  of  estimated  residual 
variance. 


These  results  show  that  deviations  from  the  actual  structure  in  crustal 
models  used  for  event  locations  contribute  significantly  to  location  errors. 
We  therefore  recommend  that  future  effort  be  directed  towards  improving  model 
structures  rather  than  forcing  observed  data  to  fit  some  given  model. 
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